Asymmetric Thermal Lineshape Broadening in a Gapped 3-Dimensional 
Antiferromagnet - Evidence for Strong Correlations at Finite Temperature 
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It is widely believed that magnetic excitations become increasingly incoherent as temperature is 
raised due to random collisions which limit their lifetime. This picture is based on spin- wave cal- 
culations for gapless magnets in 2 and 3 dimensions and is observed experimentally as a symmetric 
Lorentzian broadening in energy. Here, we investigate a three-dimensional dimer antiferromagnet 
and find unexpectedly that the broadening is asymmetric - indicating that far from thermal deco- 
herence, the excitations behave collectively like a strongly correlated gas. This result suggests that a 
temperature activated coherent state of quasi-particles is not confined to special cases like the highly- 
dimerized spin- 1/2 chain but is found generally in dimerized antiferromagnets of all dimensionalities 
and perhaps gapped magnets in general. 

PACS numbers: 75.10.Jm, 03.75.Kk, 03.75.Gg, 78.70.Nx 



In the conventional picture of temperature effects in 
magnetism, the excitations are long-lived at low tem- 
peratures and their lifetime decreases as temperature is 
increased. The accepted reason for this is that ther- 
mally activated excitations collide with each other lim- 
iting their lifetimes which results in a Lorentzian energy 
broadening of the lineshapes [l], Q. This model works 
well for gapless magnets with long-range magnetic order. 
Experimental investigations of the two-dimensional (2D) 
and three-dimensional (3D), gapless spin— 5/2 (S— 5/2), 
Heisenberg antiferromagnets (HAF) Rb>2MnF4 [3] and 
MnF2 [4] reveal symmetric Lorentzian lineshapes which 
broaden with temperature in excellent agreement with 
spin- wave calculations p], Q . The basic assumption be- 
hind this description is that the spin- waves interact only 
weakly and the states available to them cover an exten- 
sive region of phase space. As the population of excited 
spin- waves increases with temperature, the rate of col- 
lisions between them also increases and they fluctuate 
within this large range of states in an uncorrelated man- 
ner. The damping is then due simply to loss of coherence 
associated with the reduced lifetime of the excitations. 

The concept of thermal decoherence and the associ- 
ated Lorentzian linewidth broadening have become so en- 
trenched in current thinking that it is generally assumed 
to apply to all magnetic systems 0, [s[ . However for some 
magnets where there are strong interactions between the 
excitations and the phase space is restricted, it is not 
obvious that this reasoning should apply. Gapped anti- 
ferromagnets such as Haldane chains, spin ladders and 
dimer magnets which have a singlet ground state and 
triplet excitations are potential candidates. For the S- 
1/2, dimer HAFs a dominant antiferromagnetic exchange 
interaction Jo, pairs the spins into dimers. The ground 



state is a singlet and the excited state is a triplet gapped 
by ~Jo while interdimer interactions allow the triplets to 
become dispersive and mobile. In the case of a large gap 
compared to the excitation bandwidth, the distribution 
of dimer excitations is severely limited. In addition, the 
dimers are subject to a hard-core constraint (only one ex- 
citation is allowed per dimer), which can be modeled as 
a strong, short-range repulsive interaction. As temper- 
ature is raised and the population of thermally excited 
dimers increases we would expect them to scatter in a 
highly correlated way due to the interactions between 
them and the limited density of states available to them. 
Furthermore, the thermal development of the lineshape 
of the excitations would presumably show evidence of 
these strong correlations. 

Recently the temperature dependence of the magnetic 
excitation spectrum of copper nitrate (CN), a S— 1/2 one- 
dimensional (ID) dimer HAF (alternating chain), has 
been investigated by means of inelastic neutron scatter- 
ing (INS). This material is strongly dimerized and highly 
ID which confines the excited states to a narrow band. 
The lineshape was found to change from a symmetric 
Lorentzian to a highly asymmetric profile with increas- 
ing temperature [7]. These results were attributed to 
strong correlations between the excitations and prove 
that far from becoming decoherent with temperature CN 
behaves as a ID strongly correlated gas of hard-core 
Bosons. ID systems often present a special case in mag- 
netism where due to the geometry, collisions between the 
quasi-particles are unavoidable and interactions are often 
strong. It is important to determine whether a strongly 
correlated gas of hard-core Bosons is realizable in the 
more general case of a higher dimensional magnet where 
such collisions can be avoided. 
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Here, we investigate the temperature dependence of 
the magnons in the 3D dimer HAF S^C^Os- Although, 
the hard-core constraint also holds in this material, the 
3D nature of the interdimer interactions along with a 
much smaller ratio of gap to bandwidth results in a much 
larger available phase space than in CN. In spite of the 
fact that Sr3Cr208 represents a much more generic mag- 
netic system, we found that it also develops an asymmet- 
ric profile with increasing temperature. These results are 
important because they reveal that the conventional pic- 
ture of thermal decoherence is invalid for a much wider 
range of systems than previously believed. A strongly 
correlated gas of magnetic quasi-particles can therefore 
be realized in a whole class of systems - dimerized mag- 
nets of all dimensionalities and maybe gapped magnets 
in general. 

Sr3Cr208 consists of S-l/2 Cr 5+ ions, which form 
hexagonal bilayers that are stacked in an ABC ABC man- 
ner. The strong antiferromagnetic intrabilayer coupling, 
Jo, pairs the Cr 5+ ions into dimers giving rise to a singlet 
ground state [see Fig. Ufa)] The gapped triplet ex- 
citations are dispersive in all crystallographic directions 
due to the 3D nature of the interdimer interactions. It 
was recently discovered that S^C^Os undergoes a Jahn- 
Teller distortion that lifts the degeneracy of the singly oc- 
cupied e-orbitals of Cr 5+ and lowers the crystal symme- 
try from hexagonal (R-3m) to monoclinic (C2/c) below 
Tj T = 285 K SOS. Surprisingly THz and infrared spec- 
troscopy (Ref. 10) and Raman scattering (Ref. 11) find 
that the phonons of the monoclinic structure do not fully 
develop until much lower temperatures of ~ 120 K. Previ- 
ous single-crystal INS measurements at base temperature 
(1.6 K) revealed a gapped band of excitations dispersing 
between A = 3.4(3) meV and E upper = 7.10(5) meV with 



7.10(5) 



bandwidth B = E upper — A = 3.7meV [121 ]. For most 
wavevectors three triplet modes were observed arising 
from the three monoclinic twins that form below Tjt 
[see Fig. D^b)]. The data were successfully compared 
to a Random Phase Approximation (RPA) (red line in 
Fig. DJb)) which was used to extract the exchange con- 
stants listed in Fig. [TJa) [12|. Heat capacity and magne- 
tocaloric effect measurements of the Temperature-Field 
phase boundary around the lower critical field {H c \ = 
30.4 T) where the magnons first condense into the ground 
state, reveal a critical exponent of v = 0.65(2) in excel- 
lent agreement to the value of 0.67 predicted for the 3D 
Bose Einstein Condensation universality class [l3j. This 
confirms the 3D nature of the interdimer couplings and 
shows that the exchange interactions are Heisenberg with 
negligible anisotropy. 

Single crystal samples of S^C^Os were grown using 
the floating zone technique [l4] . The temperature depen- 
dence of the magnetic excitations were investigated by 
INS using two cold triple axis spectrometers, V2-FLEX 
at Helmholtz-Zentrum Berlin (HZB) and TASP at Paul- 
Scherrer Institut (PSI) [l5[. Both instruments were op- 
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FIG. 1: (Color online) (a) Low temperature (monoclinic) 
crystal structure of SrsC^Os showing the magnetic Cr 5+ ions 
only. The exchange interactions are labeled on the diagram 
and listed in the Table [T3 |. (b) Constant- wavevector scans 
along (1,0, L)h and (H,H,3)h measured on FLEX at 1.6 K 
are interpolated to obtain a colormap where color indicates 
the relative scattered intensity. The white circles give the 
extracted peaks positions and the red lines correspond to the 
RPA model for the three monoclinic twins assuming the listed 
exchange constants. The blue arrows point in the direction 
of the temperature-induced energy shift, as explained in the 
text. 



timized for good energy resolution. They were operated 
in W-configuration using a vertically-focusing Pyrolitic 
Graphite (PG) monochromator and horizontally-focusing 
PG analyzer along with a Beryllium filter. For FLEX 
the final wavevector was fixed to kf = 1.2 A -1 giving an 
energy resolution of 0.087 meV and two co-aligned sin- 
gle crystals (7.74 g) were used. One crystal (3.69 g) was 
used for the TASP experiment and the final wavevector 
was fixed to kf = 1.1 A -1 giving an energy resolution of 
0.059 meV. 

The colorplots in Fig. [2] show the energy scans per- 
formed on FLEX at (0, 0, 3)^ and (0.5, 0.5, 3)h as a func- 
tion of energy and temperature (h indicates the high tem- 
perature hexagonal structure). (0.5, 0.5, 3)h corresponds 
to the M-point where the three crystallographic twins are 
inequivalent giving rise to three modes, while (0,0, 3)^ 
corresponds to the T-point where the three twins are de- 
generate. At both wavevectors, the dispersion renormal- 
izes inwards towards the center of the band as the tem- 
perature is raised. Furthermore, the intensity of the scat- 
tering drops and the linewidth increases. The points on 
these colorplots indicate the centers of the peaks obtained 
by fitting the individual scans to Lorentzians. This func- 
tion can only model the lineshape at low temperatures 
and is increasingly inaccurate at higher temperatures. 
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The peaks become weaker with temperature and are dif- 
ficult to detect above T = 75 K. 

The lines in Fig. [2] give the energy renormalization 
predicted by the RPA dimer model using the exchange 
constants listed in Fig. [TJ The model follows the center 
of the excitations up to 40 K for (0, 0, 3) ^ and 25 K for the 
lowest and middle modes at (0.5, 0.5, 3)^. The agreement 
with the upper mode is poor. The inwards renormaliza- 
tion of the dispersion occurs because as temperature is in- 
creased the dimers become thermally excited which inter- 
feres with the intersite interactions, thus the dispersion 
tends towards the intradimer exchange energy Jq. Such 
a renormalization has been reported previously for other 
gapped dimer antiferromagnetsJ16| , such as: BaCuSi206 
[13, KCuCl 3 [18] and Pr 3 Tl pjfTThe strong decrease in 
integrated intensity of the mode at (0, 0, 3)h is shown in 
Fig. Efc), which also agrees well with the predictions of 
RPA. 
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FIG. 2: (Color online) Temperature-dependence of the 
constant- wavevector scans measured on FLEX at (a) (0, 0, 3)^ 
and (b) (0.5, 0.5, 3)^ compared to the RPA. The scans were 
interpolated and the color indicates the relative scattered in- 
tensity. The white circles correspond to the center of the 
peaks assuming a Lorentzian profile while the dashed ma- 
genta lines give the temperature evolution of the dispersion 
relation according to the RPA. (c) Integrated intensity of the 
excitation at (0, 0, 3)h as a function of temperature. The blue 
line gives the corresponding expectation values according to 
the RPA. 



departure from existing theory that predicts a Lorentzian 
profile which broadens symmetrically due to thermal de- 
coherence can appear asym- 
metric due to an artifact of the instrumental resolution, 
however our lowest temperature peaks are symmetric 
showing that the resolution which is temperature inde- 
pendent is good enough to prevent this effect. A simu- 
lation of the RPA dispersion convoluted with the instru- 
mental resolution calculated using the program Rescal 
(Ref. 22) confirms that a symmetric Voigt lineshape 
should be observed unless there is intrinsic asymmetry. 
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FIG. 3: (Color online) Const ant- wavevector scans at differ- 
ent temperatures measured on (a) FLEX at (0, 0, 3)^, and (b) 
TASP at (0.5, 0.5, 3) h compared to symmetric and asymmet- 
ric lineshape functions. The red lines correspond to the fit 
using Eq. [TJ where the asymmetry parameters a and f3 are 
varied. The black dashed lines correspond to the symmetric 
fit where a = f3 = 0. For the FLEX data the background 
was subtracted for clarity and the Gaussian width is fixed to 
o flex = 0.143 meV. For TASP, the Gaussian width is fixed 
to err asp = 0.079 meV. 



The individual scans at (0,0, 3)^ and (0.5, 0.5, 3)^ are 
displayed for different temperatures in Fig. [3] From these 
plots it is immediately clear that the lineshapes of the 
excitations do not broaden symmetrically according to 
standard expectations, but become increasingly asym- 
metric as temperature is raised. At (0,0, 3)^ where the 
three twins are degenerate, the mode is at 6meV which 
is close to the upper edge of the excitation bandwidth 
[Fig. [3fa)]. As temperature is increased the peak be- 
comes asymmetric with a tail extending towards lower 
energies. At (0.5, 0.5, 3)^ the lowest mode which lies 
at the bottom of the band again reveals temperature- 
induced asymmetry [Fig. [3jb)] but now the tail extends 
towards higher energies. These results are an important 



To parameterize the temperature evolution of the line- 
shape, the individual scans acquired on FLEX and TASP 
were fitted to an asymmetric peak function [7|. It is 
based on a Lorentzian but where the usual argument in 
the denominator [1 + t 2 ] is replaced by the polynomial 
[1 + (t — at 2 + fit 3 )]. The Lorentzian is then convolved 
with the experimental resolution assumed to be Gaus- 
sian: 
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In this function, a is the variance of the Gaussian, d(Q) 
is the center of the peak and A(Q) is the overall am- 
plitude. 7 is the Lorentzian width and the asymmetry is 
controlled by a(Q)t 2 and a further damping term P(Q)t 3 . 
The sign of a(Q)t 2 governs whether the peak broadens 
to higher or lower energies. 
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FIG. 4: (Color online) Lorentzian width 7 and asymmetry 
terms a and /? extracted from fits of Eq. [T] to (a) the unique 
peak at (0,0, 3)^ measured on FLEX and (b) the lowest en- 
ergy mode at (0.5, 0.5, 3) h from FLEX and TASP. The dashed 
lines are guides to the eye. (c) The inset shows the peak pro- 
files from the fit of the TASP data de-convolved from the 
resolution function. 

This function was used to fit the data. The lowest tem- 
perature peaks (1.6 K) which are symmetric were fitted 
by fixing a and j3 to zero, to give a symmetric Voigt pro- 
file. The excitations have two broadening mechanisms. 
Firstly due to the instrumental resolution and mosaic 
spread of the sample which give a Gaussian profile, and 
secondly due to the finite lifetime of the magnons which 
produce a Lorentzian profile. The resolution and mo- 
saic are temperature-independent, therefore the Gaus- 
sian width was fixed to the base temperature value. 
Above base temperature a and j3 were allowed to vary 
and the best fits are given by the solid red lines in Fig. [3] 
while the fitted parameters are displayed in Fig. 3J The 
Lorentzian width 7 increases rapidly with temperature 
while j3 increases slowly. The asymmetry term a becomes 
negative with increasing temperature for the excitation 
at (0,0, 3)^ indicating that the asymmetry is weighted 
towards lower energies, while it is positive for the low en- 
ergy mode at (0.5, 0.5, 3) ^ indicating that here the asym- 
metry is towards higher energies. From these results it is 
clear that the asymmetry tends towards the center of the 
band [see arrows in Fig. H^b)]. For comparison symmet- 
ric lineshapes were also fitted to the data by setting a 
and /3 to zero for all temperatures as shown by the black 
lines in Fig. [3l It is clear that the symmetric profile is 
completely unable to model the data at and above 20 K 
at (0,0, 3)^, and at and above 15 K at (0.5, 0.5, 3)^. 



In summary, we have investigated the temperature de- 
pendence of the excitations in the 3D gapped dimer HAF 
SrsCr208 using high resolution INS. The results reveal 
an inwards renormalization of the dispersion towards the 
intradimer exchange constant (Jo) with increasing tem- 
perature. In addition, the modes are damped and their 
intensity decreases as temperature is raised. Most im- 
portantly the lineshapes of the excitations become asym- 
metric with temperature being weighted towards the cen- 
ter of the band, quite unlike the symmetric Lorentzian 
lineshape expected from spin-wave theory. The RPA 
is able to model the reduction in the integrated inten- 
sity and renormalization for temperatures below 25 K. 
On the other hand the model described here predicts 
sharp excitations at all temperatures because it ignores 
correlations between the magnons, RPA can be modi- 
fied to include damping effects 18, 23, 24} . However, it 



then predicts that, the excitations would broaden sym- 
metrically with a Lorentzian lineshape for temperatures 
T < A/&b = 40 K. This contrasts strongly with our ex- 
perimental observation that the lineshape becomes asym- 
metric for temperatures above 15 K. 

Typically the lineshapes of magnetic excitations are 
believed to broaden symmetrically with a Lorentzian 
profile as temperature increases. The mechanism is as- 
sumed to be a reduction in lifetime due to thermal de- 
coherence as the excitations scattering off each other 
@, 0, [H, S The asymmetry observed in Sr 3 Cr 2 8 
is therefore unprecedented because it suggests that far 
from becoming increasingly incoherent, the magnons in- 
teract strongly with each other and scatter in a coher- 
ent manner. The dimer excitations interact because of 
the interdimer exchange constants and the hard-core con- 
straint which can be viewed as a strong short-range re- 
pulsion. When a magnon is created at finite tempera- 
tures by the neutron scattering process, it will interact 
and scatter from the thermally excited magnons. The 
energy required to create it at finite temperature is then 
different from that required at zero temperature because 
the energy of interaction must also be considered. As a 
result the observed dispersion relation and lineshape are 
broadened predominantly in the direction of the highest 
density of states. Rather than creating individual dimer 
excitations, their correlations within the gas of thermally 
excited dimers must also be taken into account. 

There is to our knowledge only one conclusive prior 
observation of asymmetric thermal lineshape broadening 
in the alternating chain compound CN [7] . This material 
is highly ID and strongly dimer ized with a very large 
ratio of gap (A = 0.35 meV) to excitation bandwidth 
(A/B = 3.5). The asymmetry is attributed to the hard- 
core constraint and the severely restricted phase space 
due to the large gap and one-dimensionality. A theoret- 
ical model was developed applicable to this special case 
which agreed well with the data @, l25l - [27| . The asymme- 
try was also parameterized as a function of temperature 
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using Eq. [T] as described in Ref. and can be compared 
to Sr 3 Cr 2 08. After rescaling the Sr 3 Cr 2 08 data by the 
CN intradimer exchange constant and refitting, the val- 
ues for the Lorentzian width 7 and asymmetry a were 
found to be similar (within 30%) to those of CN. 

The dimers in S^C^Os are coupled in 3D and the gap 
is much smaller than for CN. The ratio of gap to band- 
width is 0.9 in S^C^Os compared to 3.5 for CN, there- 
fore, the phase space available for scattering should be 
much larger and the interactions weaker. Therefore, the 
phase space available for scattering should be much larger 
and the interactions between magnons weaker. In spite 
of this, the asymmetry of the lineshape is still clearly 
evident and of similar magnitude. 

These important results contradict the long-standing 
view that excitations become increasingly incoherent 
and short-lived as temperature increases. Instead, 
they suggest that a strongly correlated gas of magnetic 
excitations is not confined to the special case of highly- 
dimerized, ID antiferromagnets, but is realized much 
more generally in dimerized HAFs of all dimensionalities 
and gap sizes and perhaps in gapped magnets generally. 

We thank F. H. L. Essler, A. J. A. James, D. A. Ten- 
nant, F. Groitl, K. Habicht and A. M. Tsvelik for helpful 
discussions. This work is partially based on experiments 
performed at the Swiss spallation neutron source SINQ, 
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